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We investigate the implications of direct CP violation (CPV) on the determination of the unitarity 
triangle angle 7 from B — > DK decays. We show that 7 can still be extracted even with the inclusion 
of direct CPV in charm if (i) at least one of the D decays has negligible CP violation; and (ii) data 
from charm factory at threshold are used. If approximate expressions without including direct CP 
violation in charm are used, this can result in a shift in 7 that is O(ro /tb)- It is modest for 
B — > DK but can be 0(1) for B — > Dir. We illustrate the size of the shift on an example of 
Gronau-London-Wyler method. 
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INTRODUCTION 



Both LHCb pQ and CDF [2] have measured nonzero CP 
violation in the charm sector giving a combined value for 
the difference of CP asymmetries 
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\a C p = a C p{K K + ) - a C p(TT vr+) 
= (-6.45 ±1.80) x 1CT 3 . 



(1) 
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Here a/ = acp(f) is the time integrated CP asymmetry 
for D decaying to a CP eigenstate / 



a f 



_ T(D° -> /) - r(L>° -> /) 

r(D^/) + r(5^/)' 



In general af is a sum of three contributions 



dir 1 7; 

a f + a. 
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where a^ u is CP violation in the decay (or direct CP vio- 
lation), a™ CP violation from mixing and a l f CP violation 
from interference of decay and mixing [3]. Mixing effects 
are universal, independent of the final state, and are fur- 



(f)oK amplitudes, where / is a common final state to both 
D° and D° decays [HHS] (see also reviews in pHHH]). Due to 
the interference one can probe the relative phase of the two 
amplitudes, which is related to 7. In the originally proposed 
methods an important assumption was that there is no CP 
violation in D decays^] Below we relax this assumption 
and show that the effect of CP violation in charm can be 
included. For a completely general treatment additional 
information is needed. If this information is not included 
we show how big of an error is introduced when assuming 
no CPV in charm decays. 



II. GENERAL CONSEQUENCES OF CP 
VIOLATION IN CHARM 



We first focus on singly Cabibbo suppressed decays. The 
D° (D°) decay amplitudes Af (Af) to CP eigenstate / can 
in general be written as 



A f = A(D° -> /) = A } 
A f ^A(D a ^f)=Aje 



[l + r f 
T f [l + r/e^-'W] 



(0) 
(7) 



thermore suppressed by the D° — D° mixing parameters where A^e^^t is the dominant singly-Cabibbo suppressed 
x, y ~ O(10~ 2 ) and are in general negligible (the exact size 
does depend on the time interval the experiments integrate 
over, however). For instance, LHCb quotes [1] 



o-k+k- 



2 K 
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(0.10±0.01)a ind , (4) 



so that ai n d = a m + a 1 can be safely neglected. 

In the present manuscript we are interested in the effect 
that nonzero a^ lr can have on the methods for determining 
the CKM unitarity triangle angle, 7, [5] 



r 

tree amplitude. In general it can have a weak phase (fj , 
however, in the SM this is zero. Since we are primarily 
interested in the extraction of 7 assuming the SM, we will 
set (jPf — in the following. The ratio Tf denotes the rela- 
tive magnitude of the subleading amplitude due to penguin 
diagrams, while bf and <f>f are the strong and weak phase 
differences, respectively. Direct CP asymmetry is then 



dir 

a f 



2r/ sin Sf sin (ftf 
1 + 2rf cos 5 f cos (/>/+'■ 



-2rf sin 5f sin ipf, 



7 



Vu d v cb v: b v: d 

v*,v*v ub v cd 



(66 ± 12)°, 



(5) 



from B — > D^K^*' decays. These methods rely on the 
interference of B ->• D°K ->■ (f) D K and B -» D°K 



(8) 

where the last expression is valid for Tf 1 which is a 
good approximation in D decays. 

In order to have CP violation in the SM all three genera- 
tions need to participate. This means that in D decays the 



1 In 6 9 D°-D° mixing was also neglected. The effects of xjj ^ 
* mcm293@cornell.edu in the measurement of 7 are discussed in [T2HT5] . D°-D° mixing 

T jure.zupanSScern.ch will be also neglected throughout the present manuscript. 



CP violation is suppressed by the CKM elements between 
the third and the first two generations, giving a naive esti- 
mate 77 ~ 0([VcbVub/V cs V us ]a s /w) ~ 1(T 4 . The penguin 
contraction matrix elements seem to be enhanced in the 
SM, giving an estimate 77 ~ and the observed size 

of Aa C p [23 Hlj (see also [HH20]). The enhanced value 
or 77 could also be due to NP [21]. In this paper we are 
interested solely on the effect of the observed Aacp on the 
extraction of 7 from B —> DK. We can set aside the ori- 
gin of Aacp, as our conclusion will not change, as long as 
the tree contributions are SM-like and do not carry weak 
phase. 

The B~ — > DK" amplitudes are 



A(B~ 
A{B- 



D°K-)=A B , 
D°K-)=A B r B e l{5B -^, 



(9) 
(10) 



where r B = 0.099 ± 0.008, while S B = (110 ± 15)° is the 
strong phase [5 Q The amplitude for B~ — > J'dK~ decay 
is then 



A(B- —>f]jK~) = A B A T f [l + r B 



(11) 



The amplitude for the CP conjugated process is obtained 
by flipping the signs of all the weak phases, 7 — ¥ —7, <fif —¥ 



A(B+ ^f D K+) = A B A T f [l + r B e l( - SB+ ^ 



(12) 



We first discuss the errors introduced by negelecting the 
effect of CPV in charm. If 77 = 0, we have 

A(B~ ^f D K~) = A B A T S [1 + t b S Sb -^] , (13) 
A(B+ -^f D K + ) = A B A T f [1 + r B e^ B+ ^] . (14) 

Now for 77 7^ 0. Keeping just linear terms in rf,r B we 
have 

A(B~ ->f D K~) = A B A T f [1 + r B e 1 ^-^ + r/e^'-W] , 

(15) 

A(B+ — > foK + ) = A B A T f [1 + r B e l{ - SB+ ^ + r / e l( ^+^ ) ] . 

(16) 



which can be re-written as 



A(B~ ^f D K~) = A B A T f [1 + r B e^'B-^)} , (17) 
A(B+ -+f D K+)=A B A T f [\ + 



(18) 



Where, at 0(rf), 5j — (rf/r B ) cos(Sf — S B ) sin(0/ — 7). So 
the shift in 7 is 0(rf/r B ). This is the case since sensitivity 
to 7 arises at 0(r B ) therefore all corrections need to be 
compared to the size of the smaller amplitude in (13), (14). 



For each final state / there is also a shift in the S B strong 
phase, 5' B —S B = {rf/r B ) sin((5/ — S B ) cos(0/ — 7) as well as 
in the ratio of the two CKM amplitudes, which now exhibits 
direct CP violation r]^—r B = (rf/r B ) cos[<5/— 5b=f(0/— 7)]- 
In order to extract the correct value of 7 all of the above 
corrections need to be kept. Note that the CP asymmetry 
at 0(rf,r B ) takes a simple form 

A C p{B -► f D K) = 2r B sm5 B sin 7 + af, (19) 

so that the expressions can be easily corrected at this order 
by measured quantities. The CP averaged branching ratio 
on the other hand is 



Br(B -> f D K) = A 2 B Aj 2 [l + 2r B sin 6 B sin 7 



af cot(<7)], 



(20) 



so that the correction requires the knowledge of strong 
phase Sf. This phase can be measured using the charm- 
factories as we show below. For completeness we also give 
in appendix [X] the unexpanded expressions for branching 
ratios and CP asymmetries. 

Next we show that in principle one can still get 7 ex- 
actly even if there is CPV in charm decays. That is we 
can still over-constrain the system and obtain all the un- 
knowns directly from experiments. To show this we first 
count the number of parameters. Let n B be the number 
of different B decays, for instance B — > DK, B —> DK*, 
B —} Dir, and let uca and nscs be the number of Cabibbo 
allowed and singly Cabibbo suppressed D decays, respec- 
tively. Only SCS D decays are assumed to have nonzero 
CP violation, with both the branching ratios and CP asym- 
metries measured. Also, is assumed to be zero, which 
is the case in the SM. Each SCS D decay depends on four 
parameters, A J, rf,Sf, and <fif which can be reduced us- 
ing the measured branching ratio, Br(D — > /), and direct 
CP asymmetry a^p. We thus have two new independent 
parameters per SCS D decay. For each B decay we have 
three new unknowns, A B , r B , 5 B , and a common unknown, 
7. Therefore the total number of unknowns is 



Unkn : 3n B + 1 + 2nscs- 



(21) 



The total number of observables is the number of different 
B — > joK branching ratios and CP asymmetries, which is 



Obs : 2n B (n CA + n scs ). 



(22) 



2 For simplicity of notation we focus on B — > DK decay, but the 
results apply also to B — ¥ D* K and B — > DK* decays. The ratio 
of amplitudes and the strong phase difference are then rg(D*K) = 



0.121 



+0.018 
-0.019' 



8 B (D*K) = (-55±^)° andr B (DK*) = 0.118±0.045, 



5 B (DK*) = (117±^)° m 



Since this grows quadratically, while the number of un- 
knowns grows linearly, for n's high enough the system can 
be over-constrained. For instance, for n B = 2, ncA = 
^T n scs — 2 we have 12 measurables and 11 unknowns. So 
at least in principle the system is solvable. We reiterate 
that we assumed that there is no CPV in CA decays, and 
as we show next this is crucial. 
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In order to obtain general understanding of when 7 can 
be extracted let us define a CP violating phase 



a f = axg(A f /A f ). 



(23) 



To first order in rf, a/ = 2r/ cos Sf sin <f>f. For B — > DK 
decays one then has 



\A(B* -+ f D K^)\ 2 = \A B \ 2 [\A f \ 2 + r%\A f \ 2 + 

2r B \A f \\A f \cos(5 B TlTaf)]- 



(24) 



From this we see that there is an overall shift symmetry 
7 — > 7 + 0, af — > af — (j). This shift symmetry is not 
broken by measuring a* lr . Since a/ = a^ ir cot(5/, a change 
in a j can always be compensated by adjusting the unknown 
strong phase difference Sf. 

The shift symmetry has an important implication for 7 
extraction. In fact the angle 7 cannot be extracted from 
B — > DK data alone, unless we can disentangle the 7 and 
af dependence, e.g. by assuming that at least one of the 
D decays has no CP violation, af = 0. This is the case 
despite the fact that we showed above that the total num- 
ber of observables can be made bigger than the number 
of unknowns. The assumption a/ = is valid in the SM 
for the Cabibbo allowed and Cabibbo doubly suppressed D 
decays. However, the shift symmetry is so far broken only 
from the interference terms in the B decay and would thus 
require more statistics. 

One can make further progress by using the fact that 
rf <C 1 in any reasonable explanation of Aacp ■ Expanding 
to first order in 77 one then arrives at Eqs. (19), (20 1. A 



nice observation was made by LHCb collaboration in [22 , 
where they noted that a^ lr cot(Sf) appears in the same form 
for Br{B -> f D K) and Br(B -> f D n) (but with r B (B -> 
f D ir) < 1). In the ratio Br(B -> f D K)/Br{B -> f D n) 
thus the term a^ lr cot(5/) cancels. The inclusion of a^ p 
then gives only a marginal shift in 7 in the analysis of |22j . 

One can also use B — > f B K decays separately if the 
interference terms are measured using charm factories run- 
ning at "J (3770). Then the relative phases a/ — a/' can 
be measured from entangled decays. For instance, for 
#(3770) -> DD -t f D f D decays, where both D -> f and 
D — >• /' are SCS we have 



r(*(3770)^/ D /^)cxa / -a 



/' 



(25) 



where we expanded to leading power in rf,r r >. If one of 
the two final states is Cabibbo allowed (we take it to be 
/'), then 



Ai 



T(*(3770) -> f D f' D ) cx 1 - 2-^- cos(a / 



(26) 



The ratio of DCS and CA amplitudes is small, Af/Afi ~ 
A 2 ~ 0.05, while in addition the sensitivity to a/ is only 
at quadratic order, cos a/ = 1 — a 2 /2 + ■ ■ ■ . Therefore the 
^(3770) — > DD entangled decay with a CA decay on one 
side and SCS decay on the other side is not optimal for 
breaking the shift symmetry degeneracy. 



A better strategy is to use Dalitz plot decays, where 
both CA and SCS decays interfere, for instance in D° —> 
K s n + n-. In this Dalitz plot the D° -> K s p° and D° -> 
K*~ir + amplitudes interfere. We have 



T{D° -> Kstt+tt-) cx \A Ksp oBW p0 (m\,m 2 _) 



A 



K* 



-BW, 



K " 



(27) 



where BW p o K *~ are the Breit-Wigner functions, m\ = 
(Pk s +Ptt ± )' 2 i an< i the ellipses denote the other resonances 
in the Dalitz plot. For 5° decay we have similarly 

T(D° -> Kstt+tt-) cx \A Ksp0 BW p o(m 2 + ,m 2 _) 
A K , + ^-BW K ,+ (m 2 + ,m 2 _) + ■ ■ ■ | 2 . 

Because of the known strong phase variation in the Breit- 
Wigner functions one can extract the phase differences 
arg(A Ksp o/A K *- 7I +) and &rg(A Ksp a / A K * +7r - ) by measur- 
ing the corresponding interference regions in the D° and 
D° decays. Since CA decays D° -> K*-tt+ and D° -> 
K* + ir~ do not carry a weak phase (as is the case in the 
SM) we also have e,rg(A K ,+^- /A K »- 7T +) — 0. As a re- 
sult one can extract the relative phase ax sp o between the 
D° —> Ksp a and D° — > Ksp° decay amplitudes. Once 
ctK s p° is measured one can use charm-factory to measure 
the CPV phases a/ for the remaining SCS decays such as 
D -> K+K-,tt+it-,..., from *(3770) -> f D f' D decays 
with SCS D decays on both sides. 



III. GRONAU-LONDON-WYLER 
DETERMINATION OF 7 

We now show the consequences of direct CPV in charm 
on determining 7 using an explicit example - the Gronau- 
London-Wyler (GLW) method [6j[7]. We shorten the no- 
tation following [33] and define 



zi = A f A B , 
z 2 = AfA B r B e^ SB+ ^\ 
so that 

z = A(B+ - 
z' = A(B~ - 



At A 



z' 2 = A f A B r B e 



i(t>B — f) 



f D K+) = z x 
fDK~) = z[ 



Z2, 



(29) 
(30) 

(31) 
(32) 

The magnitude A B is measured from Br(B~ — > D°K~), 
A B r B from Br(B~ -> D K~), while \Af \ and \Af \ are ob- 
tained from measurements of the branching ratio and CP 
asymmetry of the D — > f decays. This fixes the magnitudes 
1^1 a| and 1 2£x,2 1'- Similarly \z\ and \z'\ are determined by 
measurements of Br(B ± — > f B K ). If there is no direct 
CPV in charm decays then \z\\ = \z[\ and \z2\ = \z' 2 \- This 
is what is assumed in the original GLW method, with the 
triangle constructions shown in Fig. [TJ The magnitudes 
1^1,21) \ z [ 2 1 5 \ z \ an d \z'\ are used to construct the two trian- 
gles (with two possible orientations), while the difference 
of the two angles, 6\ — arg(z 2 /zi) and 9 2 = a,Tg(z' 2 /z' 1 ), 
determines 7, 



27 = 9 2 - 6t = &ig(z 2 z' 1 /z 1 z / 2 ). 



(33) 
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Zl = «i 




FIG. 1. The orientation of the triangles in GLW method with CP 
violation in SCS D decays included. 



Algebraically, the difference of the two angles is given by 



2 sin 2 



1- Cl c 2 ± J(l-c 2 )(l-c|), (34) 



where we shortened c\ = cos 6\ and c 2 = cos 82 ■ 
The cosines are given directly in terms of the \zi\ from 
the two triangle relations (31) and (32). For instance 



2ziZ2\, and similarly C2 



ci = (jZir + |Z2|- 

(I4I 2 + I4I 2 -KI 2 )/|2«. 

In the derivation of ( 33 ) the assumption of no CPV in D 



decays entered at two steps 

• construction of the triangles. The angle 8± = 
arg(zi/z2) will no longer equal 5b + 7, but it will also 
contain the weak phase from charm decay amplitudes 
(and similarly for 62). 

• overlapping of the triangles. In general \zi \ 7^ \z[\ so 
that the two bases of the triangles in Fig. [I] no longer 
coincide. 

If there is CPV in charm decays, so that Af 7^ Af, then 
the relation ( 33 1 gets modified to 



i 2 - Vi = ar 



Z2Z 1 
z l z> 2/ f 

2(7- af 1 cot S f ), 



= 2(7 + a/) 



(35) 



where in the last equality we have only kept terms up to 
0(rf). The difference of the two angles, 62 — &\ is thus now 



related to 7 through a final state dependent phase shift a / 
( 23 ) . This corrects for the first of the two points above. It 



also shows explicitly that there is a shift symmetry 7 — > 
7 + </>, af — > af — <p in the problem. So unless this shift 
symmetry is broken 7 cannot be determined. 

The remaining question is how 6\ and 62 are determined. 
If they are determined from the original GLW construc- 
tion using triangles, then there is going to be an error in 
the extracted value of 7. Let's denote the phases deter- 
mined in this way as dlr and dlr . They are given explicitly 



by cf r 
with Zi 



\zi\ 2 + \z 2 \ 2 - \z\ 2 )/\2z lZ 2\ and similarly for c dlr 
z,-. The relation between the correctly and in- 
correctly determined phases is obtained after some algebra 
to be 



sin 



/A0 dir 



A0 

~2~ 



-C(0i,0 2 



x a 



f 



(36) 



where we again kept only terms up to 0(ry), and for short- 
ness of notation defined A0 dir = 6> dir -(9 dir and Ad = 2 -0i. 
The multiplicative factor in the second term is 

1 M\ - s 1 



C 01,02 



-2 



e 01,02 



4r B 



C(6 



(37) 
with 

C(0i, 2 ) = (ci - c 2 ) [l + v /(l- c 2)(l- c 2)/( ClC2 )] . (38) 

To the order we are working this function can be evaluated 
using cf" instead of Cj. 

Using GLW without taking into account CPV in charm 
one would conclude that A0 dlr equals 27, a conclusion that 
is incorrect up to the shifts shown in (35) and (36). To 
stress this let us define 



„dir 



dir 



7 -" s A6 )d.y z ( 3Q ) 

The relation between the incorrectly determined and the 
correct values for 7 are then from ( |35| ) and ( 36 ) given by 

cot (5/ - „ ; „ /o . A C(0i,02) 



dir 



7 ~ 7 + a 



,dir 
/ 



4rs sin(27) 



(40) 



or numerically, using Tb = 0.1, 7 = 68°, 



dir 



-,dir 



7 = 7 +a / 



cot 5/ - 3.6 C(0i,0 2 



(41) 



We see, that the effect vanishes in the limit of no direct 
CP violation, a dlr = 0. Since the sensitivity to 7 is propor- 
tional to t\b, part of the shift 7 — 7 dlr is relatively enhanced 
by 1/rs- This results in a shift that is a factor of a few 
times the value of a dlr for B — > DK decays. So until the 
precision on 7 does not reach a level of a few percent this 
shift could even be ignored. This is not true for 7 extrac- 
tion from B — > Dir. There the ratio of the two amplitudes 
is in this case 



r B (Dir) 



tb 



\ 2 r B - 0(0.5%). 



(42) 



V us V c d 

V cs V u d 

Therefore the corrections due to CP violation in charm on 
the 7 extraction in the case of B —> Dir decays is 0(1), 
since it is enhanced by an extra factor of 1/A 2 ~ 25, and 
needs to be taken into account. 
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IV. CONCLUSIONS 

We have provided a strategy that can be used to cor- 
rect the extraction of 7 from B — > DK and B — > Dtt 
decays for the effects of direct CP violation in charm de- 
cays. The extraction of 7 requires that there is no CP 
violation in Cabibbo allowed and Cabibbo suppressed D 
decays, while the observed CP violation in singly Cabibbo 
allowed D decays can be included using our formulas. For 
this the knowledge of direct CP asymmetry in D — > / de- 
cays, a,f, is needed, along with the knowledge of the rela- 
tive strong phases 8 / between the interfering amplitudes in 
the D decays. This strong phase can be obtained from high 
precision charm factory running at 4'(3770) in conjunction 
with CP violation measurements in D Dalitz plot analyses, 
e.g. in D — > Kstt + it~. Alternatively, it can in principle 
be determined also in B — > D^-*'K^*' decays as well, if at 
least two different types of B decays with several different 
D — > f decays are considered. Another approach, valid to 
first order in the ratio of penguin and tree matrix elements 
in D —> f SCS decays, 77, is also possible [Hj. The shifts 
due to a d j P ^ in BR{B f D K) and BR{B -S- f D n) 
are the same at C(r/). The 7 extraction can thus be ap- 
propriately modified to this order to include a^ P 7^ 0, if 
B — » JdK and B — > fpir decays are used simultaneously . 

If CP violation in D decays is ignored in the extraction 
of 7, the resulting shift is of the order of 0(a'j lr /rs) ~ few 
degrees. With increasing precision of 7 determination the 
inclusion of CPV in D decays will therefore soon become 
important. For the future 7 extraction from B Dn 
decays the shift in 7 is 0(1) so that the inclusion of the 
CP violating effects will be essential from the start. 



Note added: While this paper was being finalized Ref. 
|24| appeared, which stresses the importance of including 

dir 
a f . 
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Appendix A: Complete expressions 

In this appendix we collect the complete expressions for 
CP asymmetries and CP averaged branching ratios includ- 
ing direct CP violation in charm. The CP averaged branch- 
ing ratios are 



Br(B -> f D K) 
1 



(\A(B~ -> f D K~ 



=A%(A 



\A(B+ 
2 2 

r f r B~ 



IdK 



r 2 f + r 2 t r 2 R + 



' 1 B ' f 

+ 2t\b cos 5b cos(7 + 20j ) + 2r/ cos (5/ cos </>/ 
+ ArfTs cos<5{, cos Sf cos(7 + 0/ + 2</>J ) 
-I- 2rjr B cos 5f cos </>/ 
+ 2r^r B cos(5 B cos(7 + 20 / + 20j) , 

while the CP asymmetry is 



(Al) 



A C p(B ->■ f D K) 



A%{A T f f 
By(B -> f D K) 



2tb sin 8b sin(7 + 2<^)J) + 2r/ sin Sf sin 0/ + 4r/rs sin 5b cos (5/ sin(7 + <f>f + 2</>J) 
2rfr B sin Sf sincjjf + 2r 2 rs sin Sb sin(7 + 2cj)f + 2</>J) , 

(A2) 



where Br(B — » JdK) is given in ( Al ) 
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